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We present our analysis of B physics quantities using non-perturbatively matched Heavy Quark 
Effective Theory (HQET) in = 2 lattice QCD on the CLS ensembles. Using all-to-all propa- 
gators, HYP-smeared static quarks, and the Generalized Eigenvalue Problem (GEVP) approach 
with a conservative plateau selection procedure, we are able to systematically control all sources 
of error. With significantly increased statistics compared to last year, our preliminary results 
are m b (m b ) = 4.22(10)(4),GeV for the MS b-quark mass, and / B = 193(9) stat (4) 2 MeV and 
/b s = 219(12) sta tMeV for the B-meson decay constants. 
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Figure 1: An illustration of the observed tension between the determinations of |V u b | from leptonic, exclusive 
semileptonic, and inclusive semileptonic channels; ±1(7 bands are shown. 



1. Introduction 

Weak decays of heavy mesons constrain the CKM matrix encoding the flavour-changing weak 
interactions. Besides experimental data, lattice QCD results for low-energy hadronic matrix ele- 
ments decisively contribute to precision tests in the beauty sector. Since the significance of these 
tests is predominantly limited by theoretical uncertainties, lattice computations with an overall 
accuracy of a few percent are highly desirable. 

A quantity that has attracted much attention in this context recently is the CKM matrix element 
|V U b|. It can be determined in different ways, in particular from inclusive semi-leptonic processes 
B — > X u £v (whose theoretical treatment involves the use of the optical theorem, the heavy quark 
expansion and perturbation theory), from exclusive semi-leptonic decays B — > iziv involving the 
hadronic form factor f+{q 2 ), and from exclusive leptonic decays B — > TV involving the hadronic 
decay constant /b (where the latter two cases require lattice determinations of f+(q 2 ) and /b hi 
order to extract | V u b|)- Currently, there is a ~ 3a tension between the two exclusive (semi-leptonic 
and leptonic) determinations of |V" u b|> as well as tension with the estimate from inclusive B -meson 
decays. Precise and reliable lattice calculations with good control of systematic errors are required 
in order to answer the question whether this tension really hints at New Physics in the B-sector. 

In order to describe B physics on the lattice with controlled errors, it is necessary to cover 
multiple physical scales differing by several orders of magnitude; in principle it would be de- 
sirable to have Air = L _1 <C m K , . . . , <C cT x = Atjv- What is achievable in practice is 
L > 4/nijz 6fm to suppress finite-size effects for light quarks, and a < 1/(2otd) ~ 0.05 fm to 
tame discretization errors in the charm sector and allow the use of relativic c-quarks. However, 
the b-quark scale ~ 4ra c has to be separated from the others in a theoretically sound way be- 
fore simulating the theory; in this work, we employ the Heavy Quark Effective Theory (HQET) 
formulation [1] for the b-quark in heavy-light systems. 

In this update on our earlier analysis [2], we are now able to include more statistic, as well 
as new measurements including a strange valence quark (although not all of the latter are fully 
analysed yet). The newest analysis also benefits from an improved estimate of the matching scale 
L\ = 0.401(13)fm, which previously dominated our eror on rat,, as well as a refinement of the 
procedure used to select the plateau regions from which we extract our estimates of the large- 
volume observables. 

2. Methods 

2.1 Non-perturbative HQET at 0(1 /ra b ) 

By performing a systematic asymptotic expansion of the QCD Lagrangian in AQco/'Wb *C 1, 
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and truncating 0(Ag CD /m^) contributions, one arrives at the (continuum) HQET Lagrangian: 

-^HQET (x) = Wh (*) D Wh (x) - OJkin fim (x) ~ ft) sp in ^spin (x) , (2.1) 
^kin (X) = Wh (X) D 2 Yh (x) , ^spin (x) = y h (*) (7 • B ^ h (x) . 

Similarly, the time component of the heavy-light axial current Aq (at p = 0) expands as 

,s . s 

, Ajf 1 = W\ 7o7s Vh , A { 1] = v?] 7 5 7; Vk- (2.2) 

In the HQET approach, the 1/m-terms appear as local operator insertions in correlation func- 
tions, and therefore the renormalizability of the static theory carries over to HQET (at each or- 
der in \/m). This ensures the existence of the continuum limit, once the HQET parameters 
ft), G {wibare ' Z^^ 7 > c a ■> ' ^spin} have been fixed through non-perturbative matching [3] so that 
no uncancelled power divergences in a~ l (which are induced by operator mixing in the effective 
theory) remain that would spoil taking the continuum limit. 

The strategy of non-perturbative matching [4, 5] is as follows: The matching is performed 
in the Schrodinger Functional scheme in a small volume L\ s» 0.4 fm, where due to arm, *C 1 
simulations with a relativistic b-quark are feasible. The HQET parameters ft), are fixed by imposing 
the matching conditions 

<J>HQET (z a) = £QCD ( Q) d,QCD (Z)()) = limO.QCD^^) (2 3) 

so that ft), inherit the quark mass dependence from non-perturbatively renormalized QCD via their 
dependence on z = L\M, where M is the RGI quark mass [6]. We then use a recursive finite-size 
scaling procedure to perform the step L\ — > L2 = 2L\ and finally to make contact to physically large 
volumes L«, > max(2 fm,4/m^). 

As a result of this procedure [5], the Nf = 2 HQET parameters (Oi(z,a) (which absorb the 
power divergences of HQET) are non-perturbatively known for a number of z-values around the b 
mass at the lattice spacings used in our large-volume simulations. 

2.2 Large volume computations and techniques 

Our large-volume measurements are performed on the Nf = 2 CLS ensembles, which use the 
plaquette gauge action and non-perturbatively 0(a) improved Wilson quarks and were generated 
using the DD-HMC [7] and/or the MP-HMC [8] algorithms. The ensembles fulfill the condition 
Lm n > 4, and we use a range of pion masses (190 < m K < 440) MeV at three lattice spacings 
(0.05 < a < 0.08) fm, where the scale has been set through /k [9]. 

In the computation of the static-light correlation functions, we use HYP smearing for the static 
quarks [13] as well as a variant of the stochastic all-to-all propagator method for the relativistic 
quarks (with multiple noise sources per configuration and full time-dilution) [10, 12] in order to 
improve statistical precision. 

To control excited state contaminations to the HQET energies and matrix elements, we solve 
the Generalized Eigenvalue Problem (GEVP) [11, 12] 

C(t)Vn(t,t )=X n (t,to)C(to)Vn(t,to), t < t < 2t Q , (2.4) 



.HQET _ 7 HQET 
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a [fm] 


7" 3 v T 

L, A 1 


111% 1 IVlc V 1 


it 


5.2 


0.075 


32 3 x 64 


380 


1000 






32 3 x 64 


330 


500 


5.3 


0.065 


32 3 x 64 


440 


1000 






48 3 x 96 


310 


500 






48 3 x 96 


270 


600 






64 3 x 128 


190 


600 


5.5 


0.048 


48 3 x 96 


440 


400 






48 3 x 96 


340 


900 






64 3 x 128 


270 


900 



Table 1: Overview of the CLS configurations used in this analysis. The ensembles shown in bold are new 
to the present analysis. 



for an N x N correlator matrix C(t) with N = 3, and derive the energies and matrix elements 
from the eigenvalues and eigenvectors X n , v„, such that the corrections to the energies behave like 
c< exp{ — (£at+i —E\)t} and the corrections to the matrix elements like oc exp{— (E^+i — E\)to}x 
exp { — (£2 — E\ ) (t — to) } ; for details see [11,2]. We minimize our systematic errors by a conserva- 
tive choice of plateau ranges: for fixed ? max we vary t m i n such that the errors o(t m i n ) of the plateau 
average A(f m ; n ) fulfil a sta t > 3a sys . Specifically, for each value of ? m j n we compute 

/, \ _ \A(twin)-A(twm-8)\ 

'{'mm) — / — ~- r =-, jtt- ! 

( t mm) + 2 {t min -8) 

where 8 = |ro ~ 2/(En + \ — E\ ) is chosen such that we expect the influence of excited-state con- 
tributions to have decayed by a factor ~ e 2 , and take the first value of f m j n satisfying r(? m ; n ) < 3. 

3. Results 

Combining the HQET parameters with the GEVP results for matrix elements and energies, we 
obtain observables depending on the pseudoscalar (sea) mass m n and lattice spacing a, as well as 
(through the HQET parameters) on the heavy quark mass parameter z. 

3.1 Mass of the b-quark 

We fix nib by imposing m^{zb,in^ v ,a = 0) = m^ p = 5279. 5MeV through the fit ansatz 

m B {z,m JC ,a,UYPn)=B(z) + Cml--^-,ml+DuYPna 2 , g = 0.51(2) [16]. (3.1) 

Using the NLO mass definition of HQET, m B = m bare + £ stat + (0^^° + co spin E spin , we find 

Z b = 13.34(33)(13) z , or equivalently m b (m b ) = 4.22(10)(4) :; GeV . (3.2) 
Having fixed the physical mass of the b-quark, we interpolate the HQET parameters to z = Zb- 
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Figure 2: Examples of preliminary results from an ensemble with a = 0.048 fm, ~ 340 Me V and L 3 x 
T = 48 3 x 96; left: HYP2 static-light; right: HYP2 static-strange. Shown are the effective masses or matrix 
elements as a function of t = to + a, together with a band indicating the extracted plateaux. The result of a 
global fit to the data at to /a > 6, which also includes a large number of data points not shown here (including 
results from t > to +a, and from N = 4,5), is also displayed as a solid curve, showing good agreement 
between the fit and the more conservative analysis ultimately employed. 



3.2 B-meson decay constants at NLO of HQET 

We determine the B r -meson decay constant for light (r = d) and strange (r = s) quarks through 

ln(« 3 / 2 / Br = ln(Z^ QET ) + H^PT) + + CO^p^ + O^pf* + ^ pf . 

(3-3) 

In order to estimate a systematic error in our combined chiral and continuum extrapolation, we use 
both a fit motivated by Heavy Meson Chiral Perturbation Theory (HM^PT) [14, 15] and a linear 
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Figure 3: Left: HMjPT extrapolation of /b; centre/right: linear extrapolation of /b and /b s . The blue, 
red and green points correspond to ensembles at a = 0.075 fm, 0.065 fm and 0.048 fm, respectively. Filled 
symbols denote HYP2, empty symbols HYP1; the joint continuum and chiral extrapolation is shown in 
black, with the fit formulae evaluated at each given lattice spacing shown in colour (solid for HYP2, dashed 
for HYP1). 



fit in m|: 



M (% , a, HYPn) = b r + c r m 2 + <i r , H YPn a 1 , 

3 l + 3g 2 



f B (m K ,a,UYPn)=b' 



1 



4 (4nf K ) 



m n \n(m n ) 



(linear) (3.4) 
+ c'm\ + 4 YPn a 2 . (HM*PT) (3.5) 



Both of these fit formulae are based on a simultaneous expansion in a and 1/mt,, where 0(d) 
discretization effects are dropped, since they are also 0(1 /m\,); for fy, only the non-analytic terms 
coming from the lowest order in 1 //Mb are retained. 

Our analysis for Lattice 2012 (where for /b s not all ensembles are analysed yet) gives 

/ B = 193(9) stat (4) z MeV, / B> = 219(12) stat MeV. 

The quenched value was /b s = 216(5) (using ro = 0.5fm), indicating small quenching effects. Note 
that our present Nf = 2 estimate of /b is about one standard deviation larger than the previous value 
from [2], due both to increased statistics and the inclusion of additional sea quark masses. A further 
improvement is an improved control of systematic errors from excited-state contributions (which 
enlarges the statistical errors). 



4. Conclusions 



We have presented results with a significant increase in statistics and improved control of 
systematic errors. We obtain a value for /b that is similar to those obtained by other collaborations 
[17, 18, 19, 20]. In the case of /b s> we find quenching effects to be undetectable, while for /b such 
a comparison is not meaningful due to the pathological chiral behaviour of the quenched theory. 

We are now finalizing our analysis, and more detailed publications are forthcoming. We are 
also investigating further spectral quantities within our approach, such as the B-meson spin split- 
tings, and are preparing for the determination of B — > light semileptonic form factors [21]. 
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